Galactic cold dark matter as a Bose-Einstein condensate of WISPs 
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We propose here the dark matter content of galaxies as a cold bosonic fluid composed of 
Weakly Interacting Slim Particles (WISPs), represented by spin-0 axion-like particles and 
spin-l hidden bosons, thermalized in the Bosc-Einstein condensation state and bounded 
by their self-gravitational potential. We analyze two zero-momentum configurations: the 

(N- 

polar phases in which spin alignment of two neighbouring particles is anti-parallel and the 

O; 

, the mean field approximation we derive the Gross-Pitaevskii equations for both cases, and, 

supposing the dark matter to be a polytropic fluid, we describe the particles density profile 
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as Thomas-Fermi distributions characterized by the halo radii and in terms of the scattering 
lengths and mass of each particle. By comparing this model with data obtained from 42 



5— i ■ 

£J} spiral galaxies and 19 Low Surface Brightness (LSB) galaxies, we constrain the dark matter 

particle mass to the range 10 -6 — 1 eV and we find the lower bound for the scattering length 
■ to be of the order 10~ 14 fm. 



PACS numbers: 98.80.Cq; 98.80.-k; 95.35.+d 



INTRODUCTION 



^ ■ One of the biggest challenges of modern science is to determine the constituents of an unusual 

sort of matter that is at the moment only observable through its gravitational interaction with 
usual (barionic) matter. This so-called dark matter corresponds to about 27% of the total density 
energy of the universe at the present time [it], and may dominate the total mass in galaxies. 

Amongst the many proposed candidates we can point the WIMPs (Weakly Interacting Massive 
Particles) [2j] . These are particles that present a very small coupling parameter to barionic matter 
and at the same time a very high mass of the order, at least, of GeV. Such a high mass poses the 
problem of the stability of these particles in our Universe. 

Another class of dark matter particle candidates, with masses in the sub-eV range, have been 
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recently proposed as WISPs (Weakly Interacting Slim Particles) [jj, 4 ] . They may include the usual 
QCD axion, axion-like particles with a similar coupling of axions but smaller masses, and low-mass 
spin-1 bosons. 

Axions, hypothetical particles proposed by Peccei and Quinn [Ejj], are scalar helds which have 
a nonzero vacuum expectation value and keep the CP invariance of the strong interactions in 
the Lagrangian that originally possesses a U{1) invariance involving all Yukawa couplings. They 
should have a mass range of 10 -6 — 1 eV. Although the experiments have failed so far to prove 
their existence due to the effect of an almost collisionless scattering with baryonic particles [(J, 
they are an interesting possibility for cosmology, since low mass axions are predicted to have been 
formed shortly after the Big Bang, and may constitute the dark matter component of the present 
Universe [?J| • Experimental searches are in execution or in the planning phase P, [^j . 

Since axions (axion-like particles) are defined as spin-zero bosonic particles, several authors 
1^ ] have suggested a cold dark matter fluid composed of a self-gravitating Bose-Einstein condensate 
(BEC) of axions (axion-like particles) constituting the galactic dark matter halo. This fluid is 
made of weakly coupled self-interacting particles and presents a huge phase space density enabling 
it to suitably describe the density profile of the galaxy, as long as convenient approximations are 
assumed. Other approaches consider noninteracting particles (see, e. g., [14I fla]). with the result 
that their masses become ultralight (~ 1(T 24 eV). 

Experiments have showed that even a spinorial gas can reach the BEC phase |16|] and this fact 
motivates the introduction of a spin-1 BEC to constitute the dark matter halo. The spin-1 WISPs 
are called hidden photons or hidden bosons and their formal derivation is particularly linked to 
the string theory framework 3, J]. The proposal of these particles as the components of the dark 
matter fluid leads to the introduction of additional parameters to model the halo's density profile, 
as we shall show latter. 

Although we do not claim to identify the particle that compose the dark matter fluid, our 
treatment allows one to relate the spin-0 particles to axions (axion-like particles) and the spin-1 
particles to hidden photons as depicted in references P, 4|. It is their spins and their masses that 
are the relevant features for the method we implement here. 

The present paper is organized as follows. In section [II] we develop the theory of the spin-0 
condensate based on reference [ll| and in section IIIII we presume the spin-1 condensate in the 
ferromagnetic and polar phases and, using an appropriate statistical analysis to the fit of 42 spiral 
and 19 Low Surface Brightness (LSB) galaxies, we obtain in section IIVI the mass and scattering 
lengths of WISPs. The discussion and conclusion of our results are in section IVl 
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II. SPIN-0 PARTICLES BOSE-EINSTEIN CONDENSATE 

In this case, we treat the condensate using a scalar mean field 0, which has the following 
Lagrangian density [ttI ]: 

£ = -^-(V0*)-(V0)-^((^)(^)-(W(^)) 

2m 2i 

- <t>*u^-\\4>\\ (l) 

where m is the particle mass, U is the external potential produced by self-gravitational effect, and 
the term \\4>\ represents a 2-body point interaction with A = 2t: h 2 a/m proportional to the s-wave 
scattering length a. 

Using the expression ([1]) in the Euler-Lagrange equations, we derive the Gross-Pitaevskii equa- 
tion for the spin-0 condensate as 

h 2 

ihd t = -—V 2 (i> + U^ + 2\\^ 2 ( j). (2) 
2m 

We regard \(j)\ 2 as the particle density of the condensate, p(x), thus we normalize it to the particle 

number N (J* |</>| 2 (f 3 x = N). Assuming that the particle number is conserved in the system, we can 

■ +ir ■ s(x) 

parametrize the wave function by 0(x, t) = e~ tfl ' y/ pCx.)e l « where [i is the chemical potential 



and the S(x) is the wave function's quantum phase 111 ]. 



The Gross-Pitaevskii equation ([2]) splits in two parts, one corresponding to the imaginary part 

V • (pv) = , (3) 

with v = — , and another one corresponding to the real part of the equation 

V(2Xp + ^-\VS\ 2 + W + U) = , (4) 
2m 

where W = — ^—^^ is the quantum potential. 

For a self-gravitating Bose-Einstein condensate, the external potential U obeys the Poisson 
equation 

V 2 U = iirGpm , (5) 

where p m = mp is the condensate mass density. In the case of a static condensate, for which v = 0, 
we can have the polytropic fluid with equation of state 

p = Kp T m = Kp]^- , (6) 
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where n is the polytropic index. Making the transformation p m = p c 6 n , with p c being a constant 
and 9 a function of the dimensionless coordinate £ defined by r = [(n + l)Kp l J n 1 /AttG] 1 ^ 2 ^, the 
Euler equation for a static fluid (j4]) becomes the Lane-Endem equation: 

1 d ( t2 d0 



The condensate density profile has a uniform behavior on the center of the system, decreasing 
towards the border. Because of this behavior, the quantum potential term contribution in the 
center is smaller than the non-linear interaction term. On the other hand, on the border, the 
contribution of the quantum potential is significant. When the number of particles is large, the 
uniform region is increased and, in this condition, we can obtain an analytical solution by neglecting 
the quantum potential term in the equation (|4|). This is the Thomas- Fermi approximation, in which 
the condensate is a fluid whose density profile is limited to a region of the space. The equation 
of state has the polytropic index n = 1, K = 2irh 2 a/m' i and the Lane-Endem equation has an 
analytic solution given by 

HO - *p . (8) 

with the appropriate boundary condition 0(0) = 1, which gives p m (r = 0) = p c . Thus p c is 
recognized as the central density of the condensate. To calculate the radius R of the condensate, 
we impose the condition p m (R(£o)) = 0, which gives £o = tt and: 



In current theories for axions, these particles are assumed to have no electric charge, but they can 
have a very small mass in the range from 10~ 6 eV to 1 eV. In reference [IS] the authors obtained 
an upper limit for the self-interacting dark matter cross-section using results from X-ray, strong 
lensing, weak lensing, and optical observation of the Bullet cluster IE 0657-56. Based on this cross- 
section, Harko et al. [19] estimated the upper limit for the scattering lenght, a < 10~ 21 m (10~ 6 /m). 
Using these data in ([9]) we estimate the dark matter condensate radius to lie in the range 10 -2 pc- 
10 7 pc. This radius range encompasses the size of dark matter halos in typical galaxies, indicating 
that the axion Bose-Einstein condensate is a viable candidate to represent dark matter halos in 
galaxies. As we will see in the next sections, observational radii data constrain the mass range 
even further. 
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III. SPIN-1 PARTICLES BOSE-EINSTEIN CONDENSATE 

To derive an effective low energy Hamiltonian of a spin-1 condensate, we introduce a spinor 
field operator ^f m (r) (where m = —1, 0, 1) corresponding to a field annihilation operator for a dark 
matter particle in the spin state |1, m). The Hamiltonian operator can be written in terms of these 



field operators as follows 



20 



H= f dr^#L(r)T m (r)* m (r) + i f j drdr' £ ^(^(^^(r, r')* m (r)* m ,(r') , (10) 

m m,m' 

where T m (r) = K m (r) + U m (r) is the kinetic energy plus the external potential energy of particles 
with spin m and Vm.m/iv, r') is the interaction potential between particles with spin m and m! . 

In reference 20[j, the author observed that the interactions between particles are different for 
distinct spins. Then the system symmetries lead the two particles vector state to be unchangeable 
under the permutation of the particles. Therefore, one must have (— 1)- F+ ^ = 1 where J- is the 
total spin of the system and J is the angular momentum between the particles. In the low energy 
dynamic of the system, we only consider the s-wave scattering caracterized by two-body collisions 
with small momentum transfers and represented by a 8 function in coordinate space. Thus J = 
and T must be even and range from to 2/, where / is the spin of the particles. For a system of 
/ = 1 bosons, we have 

Kn,m'( r ' r ') = ^( r - r ') X [Ao<Wm',0 + A 2 (5 m+m ',2 + S m+m >- 2 )] , (H) 



where A^ = ^^,with ajr being the scattering length between particles of total spin T . Likewise, 
the relation Fi -F2 = (#m+m',2 + ^m+m',-2) — 25 m+m / ) o, where F is the angular momentum operator, 
and the completeness relation 1 = 5 m +m',o + ^m+m',2 + ^m+m',-2> we have the general form of this 
interaction, 

V m , m > (r , r') = 5{r - r') [c + c 2 Y l ■ F 2 ] , (12) 

where c = A ° + 3 2A2 and c 2 = ^f^. 

Finally the effective low energy second quantized Hamiltonian is 

H = j dr^*t,(r)T m (r)* m (r) + 

m 

+ |/ dr E ^ t m (r)^L(r)l' m (r)^(r) + 

m,m' 

+ |/ dr E *!n(r)*L'(r)Fi • F 2 * m (r)* m ,(r) . (13) 
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As the interatomic interaction is elastic and conserves the total number and total spin of parti- 
cles, the hamiltonian is invariant under the global U(l)^ = {e^l|</> 6 [0, 2tt]} gauge transformation, 
the S'0(3)f = {e~ lFzOL e~ lFy P e~ lFx1 } rotation in spin space where (a, /3, 7) are the Euler angles, and 
time reversal Q = {1,T} where T is the time-reversal operator. 

At low temperature, there are phases regarding to the broken symmetries of the ground state. 
To describe these phases, we break the global gauge and the rotation in spin space symmetries 
substituting the fields operator ^> m by the Bose condensate \l/ m (r) = (^ m (r)) in the Hamiltonian 
(|13p obtaining the energy functional H($ m ) = {H), and construct the Lagrangian density as 



C = -|EJft)( 9 t*m) - (d t ^* m )(^ m )] - H(* m ). Using the Euler-Lagrange equations, we 
obtain the Gross-Pitaevskii equation for the spin-1 condensate, 

h 2 v 2 



in 



dt 



2M 



+ U(r) 



^ m + c n(r)^ m + c 2 5^<F) • f mm '*m' , (14) 



where n(r) = (h) = J2m l^( r )m| 2 is the total condensate density, (F) = (F x ,F y ,F z ) is the average 
of the angular momentum components, and f mm ' = ^mF^^, is the angular momentum projection 
of state Vl/ m . For the sake of simplicity we assume equal masses M for different spin particles and 
equal external potential, U m (r) = U(r). 

Considering the total number of particles as being conserved, we can write *$> m (r,t) = 
i^ m {r)e~ % ^ t ^ 1 where Vm( r ) is a steady state and \i is the chemical potential, and therefore we 
get the time-independent Gross-Pitaevskii equation 

h 2 v 2 



2M 



+ U(r) 



i> m (r) + c n(r)^ m (r) + c 2 ^J(F) • f mm 'ip m '(*) ■ (15) 



The equation (| 15 j) had already been derived in [20] , where it was solved by the use of the 
energy balance. In the present work, we take a different approach and solve this equation by 
considering that the phases are related to the spin rotational. In order to identify these phases, 
let the condensate be rewritten by ipm{^) = \J n(r)^ m (r), where ( = (Ci? Co? C— i) T is a normalized 
spinor that transform by the spin rotation D(a,/3,j) = e~ tFza e~ lFv ^ e~ lFx,y . Kawaguchi and Ueda 



showed in [22( that there are two possible phases for the ground state related to the inert states 
that have continuous isotropy groups. There are the ferromagnetic state where £ = (1,0, 0) T and 
the polar (antiferromagnetic) state where £ = (0, 1,0) T depends on the signal of c 2 . The Physics 
of the phase diagram is simple: (F) is zero for c 2 > (i.e. A 2 > Ao) and the ground state is 
polar (antiferromagnetic), or (F) is maximal for c 2 < (i.e. A 2 < Ao) and the ground state is 
ferromagnetic. Hence, there are two distinct cases: 
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(I) Polar state, where the spinor £ and the density n(r) in the ground state are: 



p i0 

C = -= {e~ ia sin /3, cos /3, e ta sin /3) T , (16) 
v2 

n p (r) = -[ f x-U(r)-W(r)], (17) 
Co 

where M^(r) = ^37 ^-^p is the quantum potential. 

(II) Ferromagnetic state, where the spinor £ and the density n(r) in the ground state are: 

c = e ^-7) (e — cog 2 ^ ^2 CQS | sin ^ eia sin 2 ^ (lg) 

n f (v) = U^-U{v)-W(v)}. (19) 

Considering that the external potential obeys the Poisson equation and the gas is composed of 
a large number of particles, we perform the Thomas- Fermi approximation neglecting the quantum 
potential in equations (|17|) and (|19|) . Thus, we have 

V 2 n p (r) = n p (r) , (20) 

c 

V 2 n f (r) = n/(r) , (21) 

which constitute the polytropic gas with index n = 1 and the coefficient K = 2cq for polar state 
and K = 2A2 for ferromagnetic state. Using the parametrized analytic solution of Lane-Endem 
equation as density profile, we can fit the condensate radius R using the interaction parameters, 
ao and 02, and its mass. For the polar phase we have 



Rp ~ ^ 3GM3 (22) 

and for the ferromagnetic phases 

R ' = ^Gm- (23) 

If the dark matter condensate in the galaxy is in the ferromagnetic phase, i. e. Oq > a^, the halo 
radius depends only on the scattering length and the treatment is similar to the spin-0 model. 
However, if the dark matter condensate in the galaxy is in the polar phase, i. e. Oq < a 1 ^, the halo 
radius depends on the scattering lengths ag and a\. In this phase, there are two parameters to be 
determined and using astronomical data we can constrain the mass and scattering lengths of these 
particles, making use of a Maximum Likelihood analysis, as shown in the next section. 
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IV. STATISTICAL ANALYSIS 



In order to constrain the mass of the dark matter particle and its relevant scattering lengths, 
we use here a subset of the astronomical data presented in 23J]. This subset presents information 



on 42 dark matter dominated spiral galaxies (/dm > 0.5, using Chabrier initial mass function), 
studied by the weak lensing method. We are interested in these galaxies' halo radii, which range 
from 1.39 kpc to 20.09 kpc. 

Using these previously mentioned data, we can construct the Likelihood function C for the halo 
radius R(a) trough: 

N 



exp 



4=1 



1 



2 a 



[R(a) 



(24) 



where R(a) represents the theoretical radius obtained from eqs. ([9]), (122j) or f)23|) . T{ are the data 
taken from observations and cr, are the errors associated with these measurements. The errors were 
not available in the mentioned work, hence we have decided to overestimate these quantities and 
assume a worst-case scenario in which they would have half the value of each radius measurement. 

For the particle mass, for both the spin-0 and spin-1 cases, we have chosen at first to use the 
lower bound of the mass range for the axion, 10 -6 eV (see 

For the spin-0 case, in which there is only one scattering length a to be determined, the proba- 
bility density function constructed is shown in the left panel of figure (JT]). 
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FIG. 1: Probability density function for the scattering length (in fm) using 42 spiral galaxies radii and a 
dark matter particle mass of 10" 6 eV. The left panel shows the spin-0 case, and the right panel shows the 
spin-1 case. 



The peak of the probability density function constructed from the Likelihood gives us the best 
fit for the parameter in investigation. We can see in figure (pQ) that the best value is around 
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10 fm. Since we have overestimated the errors used in calculations, the width of this curve is 
also overestimated. 

For the spin-1 case, we have to make some additional assumptions. We assume that the ferro- 
magnetic phase scattering length and the polar phase scattering length have the same value, 
and that this does not differ from a obtained for the spin-0 case. Using the previously calculated 
value as an input, we find the probability function for the scattering length Oq as shown in the 
right panel of fig. ([1]). The best fit value for this quantity is around 10 -13 fm. 

We also performed the same analysis using slightly higher values for the mass. Fig. ([2]) shows 
the case for a dark matter particle with a mass of the order of 10 -5 eV. The scattering length in 
this case is found to be ~ 10 -11 (10~ 10 ) fm for the spin-0 (spin-1) case. 
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FIG. 2: Probability density function for the scattering length (in fm) using 42 spiral galaxies radii and a 
dark matter particle mass of 10~ 5 eV. The left panel shows the spin-0 case, and the right panel shows the 
spin-1 case. 



The case of a particle with a mass of the order 10 -4 eV is shown in Fig. ([3]). The scattering 
length in this case is found to be ~ 10 -8 (10~ 7 ) fm for the spin-0 (spin-1) case. 

Considering larger masses causes the scattering length to lie outside the upper bound referred 
to in section [nil 

Another sample, now containing 19 low surface brightness (LSB, dark matter dominated objects 
with foM > 0.9) galaxies radii ranging from 1.2 to 19.6 kpc, obtained from [24j], have been analysed 
using the same method. The corresponding plots are shown in Figs. HEJ 
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FIG. 3: Probability density function for the scattering length (in fm) using 42 spiral galaxies radii and a 
dark matter particle mass of 10 -4 eV. The left panel shows the spin-0 case, and the right panel shows the 
spin-1 case. 
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FIG. 4: Probability density function for the scattering length (in fm) using 19 LSB galaxies radii and a 
dark matter particle mass of 10" 6 eV. The left panel shows the spin-0 case, and the right panel shows the 
spin-1 case. 



The values of the scattering lengths obtained with the use of LSB galaxies do not differ from the 
previous results. A possible explanation for this behavior is that our expressions for the radius of 
the BEC fluid derived from Thomas-Fermi approximation depend only on the dark matter particle 
mass, and not on the astrophysical objects masses. Since the radii ranges are similar in both the 
galaxies sets we used, it is not surprising that the scattering lengths are also similar. 



V. CONCLUSIONS 



In this work, we have been able to develop a theory for a spin-1 Bose-Einstein condensate 
composed of WISPs in terms of their scattering length and their masses. We could identify two 
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FIG. 5: Probability density function for the scattering length (in fm) using 19 LSB galaxies radii and a 
dark matter particle mass of 10~ 5 eV. The left panel shows the spin-0 case, and the right panel shows the 
spin-1 case. 
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FIG. 6: Probability density function for the scattering length (in fm) using 19 LSB galaxies radii and a 
dark matter particle mass of 10~ 4 eV. The left panel shows the spin-0 case, and the right panel shows the 
spin-1 case. 



phases of this condensation state, namely, a polar phase where the total angular momentum is 
zero, and a ferromagnetic phase where the total momentum is in the maximal projection. 

Considering also spin-0 particles (related to axions or axion-like particles), it was possible to 
use this condensate to model dark matter halos in galaxies and obtain their radii in terms of the 
same parameters. 

In order to constrain the values of scattering lengths and masses for the two possible condensate 
phases, we proceeded to a statistical analysis using a set of 42 dark matter dominated spiral galaxies 
and 19 LSB galaxies radii. Even though the set used presents a relativily small number of data, we 
have been able to limit the mass of the proposed dark matter particle to the range 10 -6 — 10 e V 
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and to find the lower bound for the scattering length to be 10 -14 fm. We recall here that the upper 
bound for this quan tity have been estimated to be at the value 10 -6 fm by the use of colliding 
clusters data in 19]. 



Using this method we were not able to determine the dark matter particle's spin, because of 
the similarity between the spin-0 and ferromagnetic phase expressions for the galaxy radius as a 
function of the scattering length. Even considering the fluid to be composed of spin-1 particles, we 
were not able to distinguish the corresponding two phases by the use of this statistical analysis, 
since they are linked by the use of as an input for a^in the construction of the probability 
functions . A more comprehensive analysis involving a significantly larger set of galaxy data may 
be necessary to establish this distinction beyond doubt. Nevertheless, one of the objectives of this 
paper is to show that this method can bring important information on the characteristics of the 
dark matter fluid. 

The study of the condensate excitations and speed of sound could also allow the determination 
of the spin state of the fluid's particles. These features can be related to the ocurrence of observable 
caustics and cusps in the galactic halo phase space [25|] and will be the subject of future work. 



We would like to point out that one interesting feature of the present approach is the use of 
macroscopic quantities (e. g., astronomically determined galaxies' radii) to gather information 
on microscopic parameters (particle masses and scattering lengths) related to quantum aspects of 
the nonrelativistic fluid. This correspondence is possible because we are treating a Bose-Einstein 
condensate. Another quantum fluid where this relation can be made is the degenerated Fermi fluid 
with half integer spin particles. This system will be treated in an upcoming work. 
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